The impact of polarized extragalactic radio sources 
on the detection of CMB anisotropics in polarization 
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Abstract — Recent polarimetric surveys of extragalactic radio 
sources (ERS) at frequencies u>l GHz are reviewed. By exploit- 
ing all the most relevant data on the polarized emission of ERS 
we study the frequency dependence of polarization properties 
of ERS between 1.4 and 86 GHz. For flat-spectrum sources the 
median (mean) fractional polarization increases from 1.5% (2- 
2.5%) at 1.4GHz to 2.5-3% (3-3.5%) at v > 10 GHz. Steep- 
spectrum sources are typically more polarized, especially at high 
frequencies where Faraday depolarization is less relevant. As a 
general result, we do not find that the fractional polarization of 
ERS depends on the total flux density at high radio frequencies, 
i.e > 20 GHz. Moreover, in this frequency range, current data 
suggest a moderate increase of the fractional polarization of ERS 
with frequency. A formalism to estimate ERS number counts in 
polarization and the contribution of unresolved polarized ERS to 
angular power spectra at Cosmic Microwave Background (CMB) 
frequencies is also developed and discussed. As a first application, 
we present original predictions for the Planck satellite mission. 
Our current results show that only a dozen polarized ERS will 
be detected by the Planck Low Frequency Instrument (LFI), and 
a few tens by the High Frequency Instrument (HFI). As for 
CMB power spectra, ERS should not be a strong contaminant to 
the CMB E-mode polarization at frequencies ^>70GHz. On the 
contrary, they can become a relevant constraint for the detection 
of the cosmological B-mode polarization if the tensor-to-scalar 
ratio is <0.01. 



I. Introduction 

The radiation emitted by extragalactic radio sources (ERS) 
at high radio frequencies, mainly synchrotron radiation from 
relativistic electrons in their jets and lobes, can be highly 
polarized, with an intrinsic degree of linear polarization as 
high as ~ 70 — 75% in homogeneous sources with an unidi- 
rectional magnetic field (see, e.g., iffl. El, for comprehensive 
discussions on the subject). 

These have to be considered as maximum values, which 
could be detected only in the case of very homogeneous 
sources with highly aligned magnetic field lines, which are 
unlikely to be observed among ERS. Actually, a much lower 
degree of total linear polarization, P, is commonly observed 
in ERS at cm or mm wavelengths (see e.g., Q, JfQ), with 
only very few ERS showing a total fractional polarization, 
II = P/S, as high as —10% of the total flux density, S. 

Nevertheless, this low fractional polarization observed in 
ERS may constitute a problem for the detection of primordial 
polarization in the Cosmic Microwave Background (CMB), 
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since ERS are the dominant polarized foreground at small 
angular scales and the CMB polarized signal is only a few 
percent of CMB temperature anisotropy. Therefore, CMB po- 
larization studies need a careful determination of the polarized 
emission from foreground sources in general and of ERS in 
particular. 

CMB polarized anisotropics can be decomposed into modes 
of even-parity (E-mode) and odd-parity (B-mode). As for the 
E-modes, generated by scalar perturbations in the primordial 
universe, they were first detected by |]5] and all the following 
observations (e.g., (6), Q) have confirmed measurements 
of CMB E-modes compatible with the concordance ACDM 
model. The more elusive B-modes are generated by ten- 
sor metric perturbations, i.e. "primordial" gravitational waves 
generated during inflation, and, according to predictions of 
inflationary models, with an amplitude directly proportional 
to the energy scale at which the inflation occurred fl8], |]9]. 
A detection of these "primordial" B-modes would provide the 
first real measure of the energy scale of inflation: thus, it would 
produce a real breakthrough in modern cosmology. But the 
detection of the B-modes CMB polarization - parameterized 
by the tensor-to-scalar ratio, T/S = r - is still a great 
challenge. In fact, the Planck mission iflOl . ifTTlPI would be 
marginally able to detect tensor metric perturbations by the 
direct detection of the primordial CMB B-mode, but only in 
the case of very high values of the tensor-to-scalar ratio of 
primordial perturbations, r>0.05 — 0.1 IPT41 . Otherwise, we 
will have to rely on proposed future space experiments, like 
COrE 1131 . specifically designed to detect CMB polarization 
by virtue of a much higher sensitivity than current ones. 

The above discussion clearly illustrates the importance of 
giving the best up-to-date estimate of the average/median frac- 
tional polarization as well as of the distributions of fractional 
polarization lfl6l . 11^, at least for each one of the different ERS 
populations ( ifTTl . lfl8l . lfl9l ) which are providing relevant 
contributions to temperature anisotropies of the CMB. 

Statistical studies of polarized emission from ERS have 
shown that at —1.4-5 GHz, the fractional polarization, IT, 
increases with decreasing flux density lfl6l . ll20l . Oil . 
ll22ll . Different explanations have been proposed of this result: 
e.g., a population change at fainter flux density |20|, |22| or a 
changing fraction of radio-quiet AGN [3 |. However, the cause 
of this increase in II is still unknown. On the other hand, 
it has been shown by [16] that this increase of the degree 
of polarization with decreasing flux density is only observed, 

'Detailed discussions of the Planck Low Frequency Instrument (LFI) 
and High Frequency Instrument (HFI) expected polarization capabilities and 
calibrations are given by [12], £DD- 



at cm wavelengths, for steep-spectrum ERS, i.e sources with 
a < —0.5 if S(v) oc v a , and not for fiat-spectrum ERS 
(a > -0.5). 

Other very recent studies of the polarized emission in ERS 
have tried to analyze the dependence of the fractional polariza- 
tion with luminosity, redshift and the source environment. The 
current, still preliminary, results show no correlation between 
the fractional polarization and redshift, whereas a weak corre- 
lation is found between decreasing luminosity and increasing 
degree of polarization ll23l . According to the analysis of highly 
polarized elliptical galaxies ll24l no differences have been 
found in the source environments between low polarization 
and ultrahigh polarization sources. This result should indicate 
that a high polarization must be a result of intrinsic properties 
of ERS. 

The outline of the paper is as follows: in Sect.|II]we briefly 
summarize the main processes giving rise to linear and circular 
polarization in ERS; in Sect.|III]we present current published 
data on the polarized emission from ERS; Sect.[lV]is dedicated 
to discuss our current results on statistical properties of the 
polarized emission in ERS; in Sect.[V]we give a short review 
on the most recent cosmological evolution model for ERS; 
Sect.|VT] presents our predictions on the contributions of ERS 
polarized radiation, given by unresolved sources, to the E- and 
B- modes; finally, in Sect. lVIIl we present our conclusions. 

II. Polarized emission from ERS 

Observations as well as statistical studies of the fractional 
polarization of ERS are very interesting on their own, and 
not only because ERS constitute a major contaminant of the 
CMB polarized signal. In fact, they are an important tool in 
active galactic nuclei (AGN) research, as they give valuable 
information on the physical quantities which determine the 
characteristics of the synchrotron radiation emitted by rela- 
tivistic electrons accelerated by homogeneous and/or random 
magnetic fields in AGN Q, 0. 

On the one hand, precise measurements of the total syn- 
chrotron radiation emitted by a radio source give an estimate 
- under some assumption, i.e. equipartition between the field 
and particle energy - of the total magnetic field strength. On 
the other hand, the degree of polarization of the radio wave 
provides information on the direction of the main magnetic 
field in the source environment and also on its degree of 
ordering. The magnetic field structure can, in turn, provide 
information on the relationship between the environments of 
the ERS and their properties. 

As already noted, compact ERS typically show a degree 
of total linear polarization of a few percent of their radio 
total intensity. Therefore, magnetic fields in radio sources are 
believed to be highly inhomogeneous, or almost without order- 
ing, although the observed non-vanishing linear polarization 
gives an indication of a certain degree of ordering of the field 
G2, M- 

The precise orientation of magnetic fields lines inside jets 
and lobes is still unknown, but theoretical arguments as well 
as observational evidence show that magnetic fields are indeed 
partially ordered. As for observations, these conclusions are 



based on measurements of the orientations of the linear 
polarization, revealing coherent structures across the images. 
From the theoretical point of view, an ordered magnetic field 
is expected when shocks compress an initially random field 
(with the field B perpendicular to the jet axis) or when such 
initial fields are sheared to lie in a plane, with B parallel to 
the jet axis (23, ED, (23, ED- 

Circular polarization (CP), which is a common feature of 
quasars and BLLac objects (or simply BLLacs), commonly 
known as blazars^, is preferentially generated near synchrotron 
self-absorbed jet cores and is detected in about 30%-50% of 
these sources 021 . Il25l . In these inhomogeneous, optically 
thick, synchrotron sources, the emission of the electron pop- 
ulation at lower energies is hidden by self-absorption. These 
invisible electrons produce Faraday rotation and conversion, 
which is the most likely mechanism capable of creating the 
observed CP in blazar sources ll32l . [26 1 . In any case, the 
measured degrees of CP are generally well below the levels 
of linear polarization ll33l . and thus negligible, at least at GHz 
frequencies. 

The change in the position angle of the linearly polarized 
radiation which passes through a magneto-ionic medium, i.e. 
the Faraday rotation, can be expressed by the rotation angle 
A0[rad]=i?M[rad/m 2 ] A 2 [m 2 ] (56), (34) experienced by the 
polarization vector, where RM indicates the rotation mea- 
sure. The rotation measure, RM, is the line-of-sight integral 
i?M[rad/m 2 ] = 0.81 J njcm- 3 ]^ [fj,G]d£[pc], n e being the 
electron density and B\\ the component of the magnetic field 
along the line-of-sight ll35ll . 

If the rotation depth is the same for all the emission volumes 
of the radio source, the net result is a rotation equal to RM A 2 
of the direction of polarization, without any effect on the 
degree of polarization. To change the degree of polarization 
there must be a variation in depth, either along or transverse 
to the line-of-sight l26ll . 

Faraday rotation is commonly observed towards extragalac- 
tic radio sources. It was first observed by 11361 . who found 
an RM of —60 rad/m 2 towards the center of Centaurus A. 
More recently, observations of the rest-frame RM in AGN 
cores and jets reported values from a few hundreds to several 
thousands rad/m 2 (e.g., (33, EH, ED). On the other hand, 
ERS dominated by emission from radio lobes show lower 
values of RMs, i.e.. from a few tens to hundreds of rad/m 2 
(e.g., (33). 

III. Polarization data on ERS at cm and mm 

WAVELENGHTS 

A. Polarization surveys at 1.4 GHz 

A large-scale polarization catalogue of radio sources is 
provided by the NRAO VLA Sky Survey (NVSS) at 1.4 GHz 
iRTI . This survey covers £1 ~ 10.3 sr of the sky with 
5 > —40°. The catalogue contains the flux density S and 
the Stokes parameters Q and U of almost 2 x 10 6 discrete 
sources with S'>2.5 mJy. Extensive analyses of these data were 

2 Blazar sources are jet-dominated extragalactic objects - observed within 
a small angle of the jet axis - in which the beamed component dominates the 
observed emission 1311 . 



carried out by [20 1 and [16|. By correlating NVSS sources 
with sources from the Green Bank 4.85 GHz (GB6) catalogue 
[42 1, statistical polarization properties have been derived for 
a sub-sample of - 30,000 ERS with S 1A gh z > 100 mJy, 
divided into steep- and fiat-spectrum sources. Steep-spectrum 
sources are found to increase the degree of polarization (II) 
with decreasing flux density: in fact, the median (mean) II 
value increases from II ~ 1.1% (~ 2%) for S > 800 mJy up 
to n ~ 1.8% (~ 2.7%) for the faintest ERS of the sample, 
at 100 < S < 200 mJy. On the other hand, flat-spectrum 
sources show a lower fractional polarization (median ~ 1.3% 
and mean ~ 2%) with no significant trend with flux density 

m. 

Two recent sub-mJy surveys at 1.4 GHz with polarization 
measurements, i.e. the Dominion Radio Astrophysical Obser- 
vatory (DRAO) Planck Deep Fields project 0, ED, (23] 
and the Australia Telescope Low-Brightness Survey (ATLBS) 
1 22 1, seem to confirm the above result of a higher level of 
fractional polarization in steep-spectrum sources with fainter 
flux densities. However, there is not a clear explanation for 
the origin of the anticorrelation between II and flux density 
in this source population. 

In Fig. [T] we plot differential number counts at 1.4 GHz for 
the polarized intensity, computed from the surveys presented 
above. The continuous curve represents a fit to total number 
counts of AGNs from [44|, whereas the dotted curve is ob- 
tained from the previous fit by assuming a constant fractional 
polarization of n = 3.3% for all ERS. This value is chosen 
in order to fit NVSS data; the fit is extremely good down 
to polarized intensities of P ~ 3 mJy. At fainter P fluxes the 
predicted curve bends down whereas observational data keep a 
flatter shape, in agreement with a higher fractional polarization 
for ERS at fainter flux density levels. 

B. Multiwavelength samples of polarized ERS 

1 ) Linear polarization in a sample of the B3-VLA survey: 
Polarization measurements of sources in the B3-VLA survey 
B31 were carried out by ||46] at 10.5, 4.85 and 2.7 GHz using 
the Effelsberg 100-m telescope. Taking into account only 
sources with Sw, zghz > 80mJy, a sample of 106 objects (out 
of 208) were defined with detected polarization. They found 
that flat-spectrum sources are significantly less polarized than 
steep-spectrum ones at 10.5 GHz. The latter ERS population 
shows a median fractional polarization that strongly depends 
on the frequency, from — 2% at 1.4 GHz to — 6% at 10.5 GHz, 
indicating a high Faraday depolarization at lower frequencies. 
Moreover, they notice that compact steep-spectrum sources 
exhibit much stronger depolarization than non-compact ones, 
and that sources showing larger linear size are more polarized. 
On the other hand, flat-spectrum sources are characterized by 
almost constant and low degrees of polarization (~ 2.5%) over 
the whole wavelength range considered. 

2) The Australia Telescope 20 GHz (AT20G) survey: The 
AT20G survey is a blind survey of the whole Southern sky 
at 20 GHz with follow-up observations at 4.8 and 8.6 GHz 
earned out with the Australia Telescope Compact Array 
(ATCA) |4). The full source catalogue includes 5890 sources 




S(Jy), P(Jy) 

Fig. 1 . Differential number counts at 1 .4 GHz as a function of total intensity 
(upper points) and of polarized intensity (lower points: black empty circles 
from NVSS; red empty diamonds from DRAO Planck Deep Fields; green 
squares from ATLBS). The continuous curve represents the number counts of 
AGNs from [44], whereas the dotted curve is obtained by the previous counts 
and assuming a constant fractional polarization of 3.3% for all ERS. 



detected above a flux-density limit of 40mJy. The AT20G 
catalogue is found to be 91% complete above 100 mJy and 
79% above 50mJy in regions south of declination —15°. 
Polarization was detected at 20 GHz for 768 sources, 467 
of which also have simultaneous polarization detections at 
5 and/or 8 GHz (upper limits on the polarized flux density 
are provided for non-detections, as described in |@]). Taking 
into account upper limits, a median (mean) fractional polar- 
ization of 2.6% (2.7%) at 20 GHz is found for sources with 
S20GHZ > 250 mJy, whereas an average 17% depolarization 
is observed at 5 GHz with respect to 20 GHz l47l . Compact 
sources with detected polarization are separated into flat- 
and steep-spectrum sources: the mean values of n are 2.9% 
and 3.8% respectively, thus confirming the higher fractional 
polarization in steep-spectrum sources. 

3) VLA polarization measurements ofWMAP point sources: 
Using the VLA, ll48l carried out polarization measurements 
at 8.4, 22 and 43 GHz of a complete sample of ERS brighter 
than 1 Jy in the 5-year WMAP catalogue and with declinations 
north of —34°. The sample consists of 203 objects: polarized 
emission was detected for 123, 169 and 167 sources at 8.4, 22 
and 43 GHz respectively (at 8.4 GHz only a subset of 134 were 
observed) and 105 were detected at all the 3 frequencies. An 
accurate analysis of the statistical properties of the polarized 
intensity of the sample is done by [40|, including an analysis 
of the correlations between the fractional polarization and the 
spectral indices. Here below, we summarize the main results 
obtained by them: 

> The distribution of the fractional polarization varies 



slightly as a function of the frequency. Including sources 
undetected in polarization, the mean fractional polariza- 
tion is (II) = 2.9 3.0 and 3.5 percent at 8.4, 22 and 
43 GHz. 

> No correlation is found between the fractional polariza- 
tion at 22 GHz and the intensity spectral indices, a% 2 
and a.22 (We remind readers that the sample, selected 
at 22 GHz, is dominated by flat-spectrum sources). 

• There is a significant change in the polarization angle 
between 8.4 and 22 GHz. For 45 sources, the position 
angle satisfies the A 2 dependence, and intrinsic rotation 
measures in excess of 1000radm~ 2 were observed for a 
large number of them. 

• Polarization of 71 sources were also measured at 86 GHz 
by |j49l . The fractional polarization is typically higher at 
86 GHz than at 43 GHz. 

4) VLA polarization measurements in an ACT survey field: 
In ll50ll . VLA observations in total intensity and polarization at 
4.86, 8.46, 22.46 and 43.34 GHz are presented. The sample is 
selected from the AT20G survey and consists of 159 ERS, out 
of the almost 200 sources with flux density S20GHZ > 40 mJy 
in a field of the Atacama Comsmology Telescope (ACT) 
survey. Polarized emission for about 60 sources was detected 
at all 4 frequencies, whereas the detections are 141, 146, 89 
and 59 from low to high frequencies. Fractional polarization 
distributions are very similar at 5 and 8.5 GHz, whereas a 
trend of increasing polarization with increasing frequency is 
indicated at 22 and 43 GHz. In particular, a tail of strongly 
polarized (>10%) sources is observed at 43 GHz. Data at 22 
and 43 GHz suggest also that the polarization fraction in steep- 
spectrum sources is significant higher than in flat-spectrum 
sources. 

C. Other samples of polarized ERS at high frequencies 

• [51] presented polarization observations of 250 (out of 
258) southern sources in the complete 5-GHz 1-Jy 
sample of ll52l by using the ATCA facilities at 18.5 GHz. 
Polarized flux densities were measured for 170 sources 
(114 flat-spectrum and 56 steep-spectrum), upper limits 
were set for an additional subset of 27 sources (12 
flat-spectrum and 15 steep-spectrum), and 53 sources 
were rejected (probably extended objects). The final flat- 
spectrum sample is almost complete (80%), while only 
49% of steep-spectrum sources have reliable detections. 
In the flat-spectrum sample the median fractional polar- 
ization is ~ 2.7% and the mean ~ 2.9%, and no sources 
have n > 10%. The median fractional polarization for the 
flat-spectrum sources included in the NVSS catalogue is 
about a factor 2 lower than that at 18.5 GHz. However, a 
relevant increase in the polarization fraction at 18.5 GHz 
is noticed only for sources with II<1% at 1.4 GHz. 

• [49 1 presented a 3.5 mm polarimetric survey of ERS using 
the IRAM 30 m Telescope. Their sample consists of 145 
flat-spectrum sources with 8 > —30° and flux density 
>Uy at 86 GHz. Linear polarization is detected for 76% 
of the sample (110 objects). They found that BLLacs 
(Hmedian — 4.4%) are more strongly polarized than 



quasars (H m edian — 3.1%). This result seems to be in 
contradiction with the idea that quasars should be more 
polarized at high frequencies than BLLacs because in 
the latter sources the synchrotron self-absorbed spectrum 
is mantained up to higher frequencies. A possible ex- 
planation provided by the authors comes from the recent 
evidence that the view angle of jets in quasars is smaller 
than that in BLLacs jS3l . IBH . So, if the magnetic field 
is not homogeneous along the jet, a lower fractional 
polarization level is expected from sources better oriented 
to the line of sight (i.e., quasars). 

Moreover, for those sources with detected polarization at 
15 GHz, they found that H%qghz is larger than Hisghz 
by a median factor ps 2, and about 20% of sources have 
the Use GHz /Ri5GHz ratio larger than 4. |@9] suggest 
that this increase may be explained by a combination 
of two phenomena: the 86-GHz emission comes from a 
region with greater degree of order in the magnetic field; 
or/and, at 15 GHz, emission is still affected by Faraday 
depolarization. 

• 1 55 1 looked for polarized sources in the WMAP five-year 
data, using a new technique named filtered fusion. They 
detected polarization in 13 ERS at a confindence level 
> 99% and polarized flux density higher than 300 mJy. 

IV. Statistical properties of the polarized 
emission in ERS 

In order to provide reliable estimates of ERS contamination 
to CMB anisotropy polarization measurements, we need to 
address the following main questions about polarization prop- 
erties of ERS: 1) how the fractional polarization, n, varies 
from cm to mm wavelengths; 2) the fractional polarization 
depends on the flux density; 3) how polarization properties 
change among the different populations of ERS. In this Section 
we address these questions on the basis of the observational 
data presented and discussed in the previous Section. 

A. Flat-spectrum sources 

We investigate more deeply the polarization of ERS ob- 
served in the AT20G survey. We consider the almost-complete 
sample at 8 < —15° and, when available, we use 5- and 
8-GHz measurements to separate ERS into steep- and flat- 
spectrum sources. In Table|I] we report statistical properties 
of the fractional polarization (i.e., the mean < II >, the 
median H me d and < n 2 > x / 2 ) for flat-spectrum sources as 
a function of the flux density range and the frequency. They 
are computed using the Survival Analysis techniques and the 
Kaplan-Meyer estimator as implemented in the ASURV code 
1 56], which takes into account upper limits on the fractional 
polarization for estimating the above quoted quantities. In fact, 
when polarization is not detected, an upper limit is provided. 
About 10% of sources in the AT20G sample we use here have 
measurements only at 20 GHz (see Tabled. Although most of 
them should be flat-spectrum sources, for a more consistent 
comparison with results at 5 and 8 GHz we prefer to exclude 
them from the analysis. 



TABLE I 

Median, mean and < n 2 > J / 2 of the fractional polarization for flat-spectrum sources in the almost complete sub-sample of the 

AT20G SURVEYS AT DIFFERENT OBSERVATIONAL FREQUENCIES AND FLUX RANGES. Ntot REFERS TO THE TOTAL NUMBER OF OBJECTS IN THE 
CORRESPONDING FLUX RANGE; N mu[ti TO THE NUMBER OF OBJECTS WITH 5- AND 8-GHZ MEASUREMENTS ; Nfi TO THE NUMBER OF 
FLAT-SPECTRUM SOURCES; N det TO THE NUMBER OF FLAT-SPECTRUM SOURCES WITH POLARIZATION DETECTION. 
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> 0.5 
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[0.5, 1) 






154 


93 1.54 


2.21 


2.73 
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> 1.0 






110 


93 1.90 


2.31 


2.71 


> 0.5 






264 


186 1.71 


2.25 


2.69 


[0.5, 1) 






154 


93 1.59 


2.20 


2.68 



At 20 GHz the median and the mean fractional polarization 
do not present any significant variation between the sub- 
samples defined by S > Uy and 0.5 < S < Uy. At 
the lower frequencies, < IT > and Tl me d show a moderate 
decrease at fainter flux densities. In particular, two-sample 
tests implemented in the ASURV code (e.g., the Gehan's 
generalized Wilcoxon test and the Peto & Prentice generalized 
Wilcoxon test) yield a probability of ~ 10% and ~ 1% that 
the distributions of the fractional polarization for S > Uy 
and 0.5<S'<lJyat5 and 8 GHz respectively are drawn 
from the same parent distribution (compared to a probability 
of 30% at 20 GHz). At flux densities lower than 500 mJy, the 
high number of upper limits (>50%) makes our estimates 
unreliable. 

On the other hand, a larger fractional polarization is ob- 
served as the frequency increases, with, on average, «18% of 
depolarization at 5 GHz with respect to 20 GHz, in agreement 
with results from [47 1. The ASURV two-sample tests yield 
only a 10% probability to have the same parent distribution 
for n at 5 and 20 GHz. For a better comparison, we show in 
Fig. |2] the distributions of the fractional polarization discussed 
above as obtained by using the Kaplan-Meyer estimator. These 
distributions are well fitted by a log-normal distribution with 
Timed and < n 2 > values taken from TableU (see, e.g., the 
case at 20 GHz in Fig.0)0 

In Fig. [3] we also study the correlation of Ti^oGHz with 
H4.8GHZ an d S2oghz'- no clear correlation is found between 
the fractional polarization and the flux density at 20 GHz (the 
generalized Kendall's tau test yields a probability P ~ 60% 
of no correlation). On the other hand, as expected, we see 
a strong correlation between the fractional polarization at 20 
and 4.8 GHz (with a P < 0.01% of no correlation). Using 
the Schmitt's method from the ASURV code, we find a 
linear regression HwGHz = 1-58 + 0.46H.4.8(3£f z . This result 
seems to indicate that only sources with very low fractional 
polarization at 5 GHz (i.e., Hi 8G# Z <2%) have a significant 

3 Fig.|2] clearly shows that the log-normal distribution, when averaged in 
each n bin, is giving predictions always compatible with the observed n 
values well inside the la level, except for very few bins at high n levels 
where the statistics is very poor. We have also verified that adopting other 
distribution functions, e.g. a truncated gaussian, we are not able to reproduce 
equally well the observed distributions of n. 



increase of n at 20 GHz (see also Fig.HJ. The large offset 
term in the linear regression we find, however, could be also 
partially due to Eddington bias in the AT20G catalogue and 
to the large number of sources with upper limits in fractional 
polarization at IT>1%. 

In Fig. [5] we plot the mean and median values of the frac- 
tional polarization as obtained from the surveys presented in 
the previous section. These data allow us to cover frequencies 
from 1.4 to 86 GHz. Although there is large scatter in the data, 
we see a general increasing trend for the median fractional 
polarization: Yl me d is ~ 1.5% at 1.4GHz, around 2-2.5% in 
the range 5-10 GHz, and finally 2-3% at v > 10 GHz. The 
mean fractional polarization has a more linear increase with 
the frequency in all the samples, with < n > varying from 
2-2.5% at 1.4GHz to 3-3.5% at v > 20 GHz. 

Most of the data presented in Fig. [5] are coming from 
samples of bright sources, typically with 5>1 Jy. Samples with 
fainter sources are the B3-VLA (S c = 80mJy) and the sample 
discussed in [50|. The latter one provides polarization for 
sources with S > 40mJy. However, since at the frequencies 
20 and 43 GHz the number of polarized ERS detected in the 
sample is less than 50 percent, we have decided to consider 
only flat-spectrum sources with S > 80mJy. Spectral indices 
are estimated using flux densities at 4.86 and 8.46 GHz. 
TablelHlreports the number of detections and the corresponding 
values of the mean, median and < n 2 y 1 / 2 of the fractional 
polarization. For the two highest frequencies of the sample we 
estimate the median in two ways: firstly, we use the ASURV 
code, by taking into account the upper limits; secondly, we 
assume a fractional polarization < 1% for those sources 
without measured polarization (values indicated in brackets in 
TablelflY The large spread between the two values at 43 GHz 
is an indication that the sample, also including upper limits, 
could provide biased values of statistical properties of n at 
these frequencies. 

If we compare results from the previous two samples (i.e., 
B3-VLA and |50|) with surveys of bright sources, we cannot 
find any clear evidence of higher fractional polarization in 
faint sources. However, larger and deeper samples of data are 
clearly required to settle the question. 

In Fig. [5] we also include the < n 2 y 1 / 2 values. This quan- 
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Fig. 2. Le/f panel: distribution of the fractional polarization at 4.8, 8.6 and 20 GHz (green dotted, red dashed and black solid histograms, respectively) for 
flat-spectrum sources with S > 500 mJy in the almost-complete sub— sample of the AT20G survey. Right panel: distribution of the fractional polarization 
at 20 GHz (black solid histogram) compared to the distribution (blue dashed histogram) produced by a log-normal distribution (see also Section lV-At with 
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2.01 and < fl 2 >= 3.84 (black dotted line; see Tabled at 20 GHz. 




Fig. 3. Left panel: correlation between fractional polarizations at 4.8 and 20 GHz for flat-spectrum sources with S > 500 mJy. Green points indicate upper 
limits in IT at 20 GHz, red points at 4.8 GHz, and blue points at the both frequencies. The solid line is the linear regressions found by the ASURV code (see 
the text). Right panel: the fractional polarization at 20 GHz as a function of the flux density at the same frequency (green points indicate upper limits in n). 



tity is important in order to estimate the angular power spectra 
of the polarized signal due to undetected ERS. At v > 20 GHz, 
i.e. the most interesting frequencies for CMB data analyses, 
sources have < n 2 > 1 / 2 ^ 4%. In ll49ll optical identifications 
are provided for ERS in the [50] sample. Over a total of 145 
objects, 107 are identified as quasars and 26 as BLLacs: we 
find n med = 3.0%, 3.6% and < n 2 > 1 / 2 = 3.8%, 4.5% for 
FSRQs and BLLacs, respectively. 

B. Steep-spectrum sources 

The number of steep-spectrum sources with polarization 
measurements becomes very small at frequencies >10GHz, 
preventing any study of the correlation between fractional 
polarization and flux density (see, e.g., Fig.|4]l. In the sub- 
sample of AT20G at 6 < —15° there are 51 steep-spectrum 
sources with S > 0.3 Jy and only 25 of them have polarization 



TABLE II 

Median and mean of the fractional polarization and 

< n 2 > x / 2 from the sajina et al. (201 1) sample. for 
flat-spectrum sources we consider only objects with 
S > 80 mIy. Median values in brackets are computed by 

ASSIGNING A VALUE OF II < 1% TO SOURCES WITHOUT A DETECTION IN 
POLARIZATION. 







Flat- 


-spectrum sources 






i/[GHz] 


N S ou 


N det 


-turned 


< n > 


<n 2 >V2 


4.86 


56 


54 


2.16 


2.47 


2.97 


8.46 


56 


55 


2.25 


2.61 


3.10 


22.46 


56 


44 


3.12 (2.88) 


3.79 


4.21 


43.34 


50 


33 


4.14 (2.64) 


4.85 


5.79 


Steep-spectrum sources 


4.86 


45 


44 


2.56 


4.10 


5.79 


8.46 


45 


44 


3.48 


4.67 


5.97 


22.46 


44 


25 


6.14 (3.40) 


7.17 


8.30 


43.34 


28 


11 


6.77 


7.18 


8.03 



detected at 20 GHz. We find that the median fractional polar- 
ization for these sources is less than 2% at all the frequencies 
(see Fig. [5] lower panel). These very small values could be 
biased due to the small number of detected sources and 
to the incompleteness of the sample at faint flux densities. 
Nevertheless, from Fig. [4] we can see a general increase of the 
fractional polarization between 4.8 and 20 GHz: the ratio of II 
at 20 and 4.8 GHz is typically close to one (< 1.5) for sources 
with H^.^ghz > 2%, but become higher for sources weakly 
polarized at 4.8 GHz. In fact, the linear regression found by 
the Schmitt's method implemented in the ASURV code yields 

n 20G ff z = 0.78 + l.UIh.8GHz- 

From Fig. [5] we can observe the strong increase of the 
fractional polarization from 1 .4 GHz to 5 GHz, where Faraday 
depolarization is probably very relevant. The increase becomes 
more moderate up to 20 GHz. The large difference in the 
median values among different samples is perhaps related to 
the small samples considered and the large fraction of sources 
with upper limits in polarization. The largest sample of steep- 
spectrum sources is provided by B3-VLA (77 sources, all with 
detected polarization) and the Ricci et al. sample (71 sources, 
15 of them with upper limits): in these samples the median 
fractional polarization is ~ 5% between 5-20 GHz and the 
mean varies from 5 to 6.5%. 

C. Number counts in polarization of ERS at 20 GHz 

We provide a first estimate of number counts in polarization, 
n(P), at 20 GHz by exploiting WMAP and AT20G polariza- 
tion source catalogues. For very bright sources, we use the 
polarized source sample detected by |55l in the 5-yr WMAP 
CMB anisotropy maps (see their Table 2). We exclude from 
number counts Fornax A, Virgo A and CentaurusA because 
they are all local objects. Then, we consider the nearly- 
complete subsample of AT20G at 8 < —15° and S > 50mJy. 
This sample allows us to compute n(P) down to polarized 
fluxes of ~ lOmJy. Number counts, that are reported in 
Tablelfllksee also Fig. [6]), have been corrected for the estimated 
incompleteness of the sample (i.e., the completeness is 0.91 
for S > 100 mJy and 0.79 for 50 < S < 100 mJy H). As 
displayed by Fig. [6] number counts in total polarization, n(P), 
are almost flat at 20 GHz, at least in the flux density range in 
which they can be estimated. 

Given that an upper limit on the fractional polarization 
H up is always provided for the ERS in this sample (see 
for more details), except for 5 very bright objects with no 
information on polarization, we also give a tentative estimate 
(TablelTTTl third column) of their contribution to the number 
counts n™ p (P) between Pi and Pi + AP by means of 

< P (P) = Y1 / n ^(n)dn, (i) 

with = Pi/ 'Sk and = p *~g . The sum is done over all 
the sources without polarization detection and Sk is the flux 
density of the k-th object. T^n) is the probability function for 
the fractional polarization: we take a log-normal function with 
n med = 2.0% and < n 2 >V2 = 2.8% (see Tableffland Fig.0 
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Number counts in polarization at 20 GHz from WMAP and 
AT20G SOURCE CATALOGUES. N det indicates the number of 
OBJECTS WITH DETECTED POLARIZATION, WHILE n u P(P) IS THE 
ESTIMATED CONTRIBUTION OF AT20G OBJECTS WITH NO POLARIZATION 
INFORMATION OR WITH UPPER LIMITS IN THE FRACTIONAL 
POLARIZATION. 
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if n < n up and zero otherwise. For the 5 objects without any 
polarization information, no upper limits is considered. 

As shown by TablelTTlland Fig. [6] the contribution of sources 
undetected in polarization (displayed by empty squares) is 
generally negligible, except for bins at very faint polarization 
levels (i.e. less than 20mJy), and for bins at P>100mJy, 
due to the 5 sources without polarization information and flux 
density between 1 and lOJy. Only for the faintest bin n" p is 
larger than the uncertainty on the number counts. Therefore, 
we can conclude that number counts in polarization estimated 
from the AT20G sample are not affected by upper limits in 
the polarized flux density. 

V. COSMOLOGICAL EVOLUTION MODELS AND 
HIGH-FREQUENCY NUMBER COUNTS OF ERS 

Early evolutionary models of radio sources [58), 11591 . 
fl7l . ll60l were able to give remarkably successful fits to 
the majority of data coming from surveys at i/<10GHz, 
and down to flux densities of a few mJy. More recently, 
ifTHl and [44 1 exploited the wealth of newly available data 
on luminosity functions, multi-frequency source counts and 
redshift distributions to provide new cosmological evolution 
models of radio sources at frequencies >5GHz and <5GHz, 
respectively. These two models are based on the determination 
of the epoch-dependent luminosity functions for different 
source populations, starting from the local luminosity function 
and by adopting different luminosity evolution laws with free 
parameters fitted from observational data. 

The predictions of high-frequency number counts of ERS 
provided by the above evolution models assume a simple 
power-law spectrum for ERS. Each source population is 
characterized by an "average", fixed, spectral index, or by two 
spectral indices (at most). This "classical" modeling has to be 
considered as a first - although successful - approximation, 




Fig. 4. Fractional polarization for steep-spectrum sources at 20 GHz as a function of fractional polarization at 4.8 GHz (left panel) and of flux density at 
20 GHz (right panel). Points of different colours have the same meaning as in Fig.[5] 



but it gives rise to an increasing mismatch with observed high- 
frequency (> 30 GHz) number counts currently available. 

Indeed, the ERS energy spectra can be quite different from 
a single power-law if analyzed in large frequency intervals. 
Different mechanisms can be responsible for this: a) a spectral 
steepening due to the more rapid energy loss of high-energy 
electrons with source age, i.e. "electron ageing"; b) a transition 
from the optically thick to the optically thin regime at high 
radio frequencies; c) at different wavelengths radio emission 
can be dominated by different components characterized by 
distinct spectral behaviors ll6D . In particular, a clear steep- 
ening at mm wavelengths is theoretically expected for radio 
flat spectra of AGN core emission 11621 . 11631 . This steepening 
has been already observed in blazars (see, e.g., [64]) and has 
also been statistically suggested by recent analyses of different 
ERS samples at v > 30 GHz E2, (66), ED, ED, fl9l . 

A first attempt of taking such steepening in blazar spectra 
into account has been done by |[T9l . In this work the spec- 
tral behaviour of blazars at mm wavelengths is statistically 
described by considering the main physical mechanisms re- 
sponsible for the emission. In agreement with classical models 
of the synchrotron emission in the inner jets of blazars these 
authors interprete the high-frequency steepening observed in 
blazar spectra as caused, at least partially, by the transition 
from the optically-thick to the optically-thin regime in the 
observed spectra at mm wavelengths. Based on the published 
models of synchrotron emission from inhomogeneous, un- 
resolved, relativistic jets ll63l . [69|, [30 1, the value of the 
frequency vu at which the spectral break occurs is estimated 
as a function of the relevant physical parameters of AGNs: the 
redshift, the Doppler factor, and the linear dimension of the 
region (approximated as homogeneous and spherical) that is 
mainly responsible for the emission at the break frequency. 

Recent high frequency data 1711 . IT721 . and in particular the 
first 1.6 Planck survey [68 1 have provided new important re- 
sults on number counts and related statistics of ERS, and they 
can be used to constrain the different possible cases/models 
featuring a spectral break in the emission of AGN jets (see 



|19|). The most successful one (indicated in lT9l as C2Ex, i.e. 
the relevant emission is emitted from a more "extended" region 
in the inner jet of FSRQs) assumes different distributions of 
the break frequency for BLLacs and FSRQs. According to this 
model, most of the FSRQs should show a bend in their other- 
wise flat spectra between 10 and 100 GHz, whereas in BLLac 
spectral breaks should be typically observed at v > 100 GHz 
(implying that the observed synchrotron radiation comes from 
more compact emitting regions). This dichotomy has been 
indeed found in the Planck ERCSC 1531. 11701: almost all radio 
sources show very flat spectral indices at LFI frequencies, 
a LFift — 0.2; at HFI frequencies BLLacs keep flat spectra 
{olhfi^j— 0.5) while most of FSRQs show steeper spectra, i.e. 
oihfi < —0.5. Moreover, the same model gives a remarkably 
good fit to all the observed data on number counts and spectral 
index distributions of ERS at frequencies above 5 and up to 
220 GHz. 

For the above reasons, hereafter we adopt the number counts 
provided by the "C2Ex" case in [19|. 

A. Number counts of ERS in polarization at mm wavelengths 

To assess the contamination due to undetected ERS in CMB 
polarization maps it is necessary to know how ma ny sources 
can be found with polarized flux density P = \J Q 2 + C/ 2 
above a given flux limit Pu m . Answering this question in- 
volves the determination of source counts in polarization. 
This estimate is quite difficult to perform directly, by the 
statistical analysis of observational data, since the polarized 
signal is weak and many samples are usually defined by 
their completeness in terms only of total intensity. However, 
this problem can be overcome by using the source number 
counts in total intensity complemented by information on the 
statistical properties of the fractional polarization n, expressed 
in term of some probability function 'P(n) |fl6l . 

4 Typically, a de-biased estimator is however used for the polarized flux 
density of point sources (as, e.g., P = y/ Q 2 + U 2 — s 2 , where s is the 
uncertainty on P; see 1501 for more details). 





Battye et al. 
Ricci et al. 
Agudo et al. 



10 

frequency (GHz) 



100 




10 

frequency (GHz) 



100 



Fig. 5. Median, mean and < II 2 > x / 2 of fractional polarization for flat- (upper panels) and steep-spectrum (lower panels) sources. Black points refer to the 
almost complete sub— sample of the AT20G survey with S c = 500 mJy for flat-spectrum sources and 300 mjy for steep-spectrum sources. Red points refer to 
data presented in (40). Green points refer to the B3-VLA sample (respect to [46], we distinguish flat- and steep-sources on the basis of their spectral index 
between 1.4 and 5 GHz; values shown in the plot are obtained using the ASURV code, taking into account upper limits in the fractional polarization). Cyan 
points refer to the sample from |50 | (see text and TablelfTl empty cyan points at 43 GHz for the steep-spectrum have to be considered as upper limits). Blue 
points refer to results from 1511 and magenta points from [49 1 (solid and empty points are for FSRQs and BLLacs, respectively). Yellow points refer to the 
NVSS survey (see |16|): for steep-spectrum sources we plot results for sources with 100 < S < 200 mjy and S > 800 mJy. 



Let us discuss this point in more detail. Polarization number 
counts n(P) can be written as: 

f°° f°° rf<? 

n(P)=N V{P 1 S)dS = N V (jl jS ) — , 

JSa=P JSo=P b 

(2) 

where N is the total number of sources with S > Sq in the 
sample, V(P, S) and V(TL, S) are the probability functions of 
observing in a source of flux density S a polarized intensity 
P and a fractional polarization IT, respectively. Assuming that 
IT is independent of S, n{P) can be determined by 

r°° 

n(P)= V(n = P/S)n(S) — . (3) 

JS Q =P h 

The probability function V(H) can be constrained from the 
observed distributions of the fractional polarization. In agree- 
ment with [40] we model 'P(n) by a log-normal distribution: 

V(IL) = J— exp J - [log(n/n med )] 2 /2 ( T 2 1 , (4) 

V27TCH1I I J 

where Timed is the median of the distribution and a 2 = 
1/2 log((n 2 )/IT 2 ne(i ). These formulas are strictly valid only 
if < II < oo. However, because of the very low fractional 
polarization observed in ERS, the upper limit of II = 0.75 can 
be effectively assumed as infinite. In Fig. [2] we compared log- 
normal functions with the polarization distributions observed 
in the AT20G survey, confirming the very good fit of the model 
with data. 



A critical point for our estimates is the variation with 
frequency of the fractional polarization observed in ERS. 
The data on ERS polarization discussed in this paper clearly 
show an higher fractional polarization at 10-20 GHz than 
at few GHz in both flat- and, more prominently, in steep- 
spectrum sources. At higher frequencies, data become scarce, 
but there are still indication of a possible further increase of 
the polarization fraction from 10-20 GHz to 40-90 GHz 11491 , 
[40 1 . For flat-spectrum sources, this increase could be due to 
the combination of two different effects: 1) the polarization 
degree actually increases with the frequency; and 2) BLLacs, 
which are observed to be more polarized than quasars [49], 
become more relevant in number at higher frequencies. 

In the following predictions we take into account both 
effects. We consider two possible values for the median and 
the dispersion (n 2 ) 1 / 2 of the log-normal distribution in order 
to provide a range of estimates for number counts and power 
spectra that could take into account uncertainties in observa- 
tional data at v > 10 GHz. Tablell VI reports our estimates for a 
more "optimistic" (lower) and a more "conservative" (upper) 
case|f| The frequency dependence of the fractional polarization 
in ERS is simply modeled by means of two different sets of 
median and r.m.s. values at frequencies below and above about 
40 GHz. Although somewhat arbitrary, this choice is motivated 
by the fact that this is the highest frequency at which multi- 

5 Please note that the choice of the two adjectives, "optimistic" and 
"conservative", has been done under the point of view that polarized ERS 
do "contaminate" CMB maps in polarization and, thus, have to be removed 
from them. 



TABLE IV 

Median and dispersion of the log-normal distribution for the 
fractional polarization, ii, as a function of frequency and 
for the different radio source populations relevant at cmb 
frequencies. two cases are considered, a more "op timist ic" 
one and a more "conservative" one (see sub-sect jv-at . 
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frequency polarization samples of ERS are available. 

• At i/<40GHz we choose a lower and higher median 
and r.m.s. values of the fractional polarization based on 
the results displayed in Fig. [5] Moreover, we require that 
number counts computed using values in TablellVlbe (at 
least partially) compatible with AT20G/WMAP counts in 
polarization. As shown in Fig. [6] the more "optimistic" 
case fits quite well observational counts, especially at 
P<50 mJy. On the other hand, the more "conservative" 
case tends to overestimate number counts at lower po- 
larized fluxes, whereas it fits particularly well current 
data at high P fluxes. This is not unexpected because, 
in this latter case, we take as median and dispersion of 
the fractional polarization for blazars the values provided 
by El at 90 GHz. 

• At v > 40 GHz, we make the assumption of an increase 
of the median fractional polarization of about 20% with 
respect to the corresponding cases at t/<40GHz. This 
choice is not firmly constrained by current data sets; how- 
ever, we have been guided by the fractional polarization 
levels of ERS observed by AID at 86 GHz which fall in 
the middle between our present lower and upper cases. 

In Fig.[6]we compare the predicted number counts, n(P), of 
ERS in polarized intensity at 20 and 143 GHz. By comparing 
the two panels of the same Figure we can appreciate that n(P) 
decreases by a factor ~ 1.5 at 143 GHz at P ~ 0.1 Jy, in 
comparison with the estimated n(P) at 20 GHz. At still fainter 
polarized intensities, P ~ 0.01 Jy, this decrement increases 
due to the fact that the contribution of steep-spectrum ERS, 
which are on average more polarized, becomes less and less 
relevant at high CMB frequencies. Finally, the right panel of 
Fig. [6] indicates a possible higher relative contribution to total 
polarized counts, n(P), coming from BLLacs with respect to 
FSRQs at 143 GHz. This is a direct consequence of the present 
analysis, which is based on the observations discussed by |49|. 

The good agreement with observational data at 20 GHz, as 
well as the reliability of our predictions on ERS number counts 
in flux density, S, give us confidence in making extrapolations 
of integral number counts, N(> P), in total linear polarization 
at Planck frequencies (see TablelVll. It is expected that Planck 
LFI will be able to detect sources in polarization down 
to ~ 200 mJy at 30 GHz and ~ 300-400 mJy at 44 and 
70 GHz, whereas Planck HFI, thanks to a better resolution and 
sensitivity, should reach polarized flux limits of P « 100 mJy. 



In both cases, however, a catalogue of only a few tens of 
compact ERS should be provided by Planck in polarization. 

The above quoted detection limits in total polarization are 
estimated by the performances foreseen for the Planck LFI and 
HFI instruments and are based on pre-launch measurements 
[12 1, lfl3l of the detectors calibration and capabilities. These 
detection limits take also into account the future application 
to the Planck data - corresponding to the end of the nominal 
mission, i.e. in January 2012 - of new detection techniques, 
specifically designed for detecting compact polarized sources 
in CMB maps (see, e.g., Il73l for a recent discussion on the 
subject). These techniques have already been applied with 
success to WMAP 5-yr maps ll55l . improving on the results 
published by the WMAP team on the same data set. 

VI. Predictions on the contribution of ERS to the 

CMB E- AND B-MODES 

As for temperature fluctuations, the analysis of CMB po- 
larization measurements is usually made by the estimate of 
angular power spectra, i.e. E- and B-mode spectra [0, ID. B- 
mode polarization, that arises only from tensor perturbations, 
is expected to be extremely weak, and even for the most 
optimistic cases foreseen by inflationary models the rms signal 
is only a fraction of /iK, less than 1 per cent of the level of 
temperature anisotropies at degree scales. On the other hand, 
polarization of foregrounds (and in particular of extragalactic 
sources) is equally shared between E- and B-modes [74|. ERS 
are therefore expected to dominate the sky B-mode polariza- 
tion at sub-degree angular scales at frequencies i^<100GHz 

nm 

By using the statistical characterization of the fractional 
polarization described in the previous Sections, we are able 
to estimate the polarized angular power spectra given by 
undetected ERS in CMB anisotropy maps. First of all, we 
assume that ERS follow a Poisson distribution in the sky. The 
contribution of clustered ERS to the angular power spectrum 
of CMB temperature anisotropy is in fact small and can be 
neglected, if ERS are not subtracted down to faint flux limits 
(the signal due to clustered ERS becomes more relevant only 
at relatively low fluxes, i.e. S < lOmJy, 11761 . 11771 ). 

It is well known that an ensemble of Poisson distributed 
point sources gives rise to a flat power spectrum of temperature 
fluctuations f78l . For a sample of sources with flux density 
below some cut-off S c , the amplitude of this white noise 
spectrum is given by 

where N and n(S) are, respectively, the total and the differ- 
ential number of sources per steradian, and dB / dT is the con- 
version factor from brightness to temperature, i.e. dB/dT ss 
lO-VKVCJysr- 1 ) ( e x - Vf/{x i e x ) and x = z//56.8GHz. 

As an analogy to the expression of the angular power 
spectrum in total intensity or CMB temperature, it is possible 
to define the angular power spectrum for the Stokes parameters 
Q and U. Because point sources contribute, on average, 
equally to Q, U and to the E-, _B-mode power spectra, we 




Fig. 6. Normalized number counts at 20 (left panel) and at 143 GHz (right panel) in total intensity (upper curves and data points) and in polarized intensity 
(lower curves and data points; blue curves are for the more "conservative" case and red curves for the more "optimistic" case). Number counts are for the 
different source populations discussed in the text: solid lines represent total number counts; dotted lines are for steep-spectrum sources; long-dashed lines 
are for FSRQs; finally, dot-dashed lines are for BLLacs. In the left panel, filled squ ares rep resents our current estimates from the AT20G sub-sample (empty 
squares include the contribution from sources without polarization detection; see Sect. lIV-Cl and empty circles are from WMAP 5-years data [ 55 ] - Data points 
in the right panel are from published data 171], [72], |68|. 

TABLE V 

Expected total numbers of ERS with polarized intensity > P Hm over the full sky at Planck frequencies. For each Planck channel, 

THE MINIMUM AND MAXIMUM VALUES HERE INDICATED REFER TO OUR PREDICTIONS CALCULATED BY THE MORE "OPTIMISTIC" AND BY THE MORE 

"CONSERVATIVE" CASES PREVIOUSLY DISCUSSED (SEE TABLeUVTi. 
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W Number of polarized ERS detected in the WMAP 5-vr data at I fc> I > 5° by [55l 
at 33, 41 and 61 GHz (P lim « 300 mJy). 



assume Cei — Cbi — Cqi ~ Cue (we generally refer to 
them as polarization spectra). Following the treatment given 
by lfl6l . we have that 

Cqi = (^yW> = (^)~Wn 2 cos 2 (20)) 

= (ff) W)(IT 2 )<cos 2 (20)> 

= 1/2 (g) ' \tf)C Te , (6) 

where the Stokes parameter Q is written in terms of II and 
of the polarization angle in the chosen reference system, 
cf>. The factor 1/2 arises because of the uniform distribution 
of polarization angles. It is easy to demonstrate that the 
cross-correlation temperature-polarization spectra are null, 
e.g. C TQ e = (S 2 IIcos(20)) = 0. 

In Eq.[3] and [6] we have assumed that the fractional polar- 
ization is independent of the total intensity of the source. 
Observations for flat-spectrum sources in the flux density 
range S^lOOmJy seem to support this assumption (see sub- 



Sect. lIV-Al ). We expect that this is maintained also at fainter 
fluxes but only if FSRQs and BL Lacs are separated into two 
different populations. On the other hand, this hypothesis may 
not be true for steep-spectrum sources: a clear anticorrelation 
between II and S is observed from data at 1 .4 GHz, whereas, 
at higher frequencies, the lack of large samples of steep- 
spectrum sources does not allow to determine it. In any case, 
steep-spectrum ERS are giving a negligible contribution to 
number counts at c > 100 GHz and, thus, this lack of 
information does not affect our current predictions. 

In Fig. [7] we present the results on ERS polarization power 
spectra for the six Planck frequencies where ERS are relevant. 
The value of (II 2 ) 1 / 2 is taken according to TablellVI for the 
different radio source populations. Moreover, we consider two 
cut-offs in flux density: S c = 1, and 0.1 Jy. The former value 
is close to the completeness limit obtained by the Planck 
ERSCS for the frequency channels v < 100 GHz (68 1; the 
latter one can be seen as a (somewhat optimistic) reference 
value for the Planck high frequency channels, or a reference 
value for future experiments. 



In Fig. [7] we also plot the CMB power spectrum for the E- 
mode and for the B-mode with a tensor-to-scalar ratio r = 
0.1, 0.01 and 0.001. In this way we can have an indication of 
the level of ERS subtraction required to allow the detection of 
a gravitational wave induced primordial CMB B-mode signal. 
ERS should not be a strong constraint on detecting E-mode 
polarization or the B-modes with r>0.01. On the other hand, 
a primordial CMB B-mode signal corresponding to lower r 
values, requires the subtraction of ERS down to flux detection 
limits of ~ 100 mJy, which will not be easy, or even possible, 
with the Planck sensitivity. 

VII. Conclusions 

In this contribution we have reviewed recently available 
polarimetric surveys of extragalactic radio sources at frequen- 
cies i/>1GHz. These data point out that the typical intrinsic 
fractional polarization of ERS is around 2-5 % of the total 
flux density, S, of the source even at frequencies as high as 
20 GHz, and that in very few objects the fractional polarization 
is II>10%. This may be due to the low degree of uniformity 
of magnetic fields in the internal part of AGN jets and in 
lobes. Faraday depolarization is probably the cause of the large 
number of sources with a very low level of polarization, i.e. 
II<1%, at GHz frequencies, which also explains the strong in- 
crease of the fractional polarization observed at ^>10GHz in 
those objects. This conclusion is supported by high or extreme 
values of rotation measure RMs> 1000 radm~ 2 observed in 
some blazars (e.g., ||37l . [40 1). 

Moreover, we have studied how polarization properties of 
ERS change from cm to mm wavelengths. For flat-spectrum 
sources a weak but constant increase of fractional polarization 
is observed, with median (mean) values varying from 1.5% 
(2-2.5%) at 1.4GHz, to 2-2.5% (2.5-3%) at 5-10 GHz and 
2-3% (3-3.5%) at 10^10 GHz. Indications that fractional 
polarization in blazars could further increase above 40 GHz 
come from the recent works by iflOl . 11491 . On the other hand, 
a significantly higher fractional polarization is typically found 
in steep-spectrum sources, especially at high frequencies 
(median is 4-5% at 10-20GHz and mean between 5 and 
6.5%). However, because of incompleteness of the samples 
and of the small number of steep-spectrum sources in surveys 
at i/>10GHz, current observations could be biased by high- 
polarized objects. 

In general, we do not find any dependence of the fractional 
polarization of ERS with the flux density at high radio frequen- 
cies. However, more conclusive evidences require larger and 
deeper surveys. Nevertheless, an anticorrelation between n 
and S in blazars is expected when very faint flux densities are 
considered (see Sect. UVl i. For the flux-density ranges covered 
by available surveys, flat-spectrum sources are dominated by 
quasars (FSRQs). These objects are typically less polarized 
than BLLacs [49|, which become increasingly relevant at 
fainter flux densities and become the dominant population at 
5<10mJy flS), ED. 

We also discuss a formalism to estimate ERS number counts 
in polarization and to predict the contribution of unresolved 
ERS to angular power spectra at CMB frequencies. As a first 



application, we attempt to predict how many polarized ERS the 
Planck Satellite will be able to detect in the different channels 
sensible to polarization measurements: we expect that only a 
dozen polarized ERS could be detected by Planck LFI, and 
a few tens at the HFI frequencies. Although the number of 
Planck detected sources is low, these data will allow us to 
study the frequency dependence of the fractional polarization 
in a wide range of frequencies, from 30 to 353 GHz, thus 
providing original and valuable information on polarization 
properties of ERS in the innermost regions of AGN jets. 

Finally, our results on polarization power spectra demon- 
strate that ERS should not be a strong contaminant to the 
CMB E-mode polarization when observing at frequencies 
^> 70 GHz. Moreover, it seems unlikely that ERS will have 
a significant impact on our ability to detect the B-mode 
polarization from primordial graviational waves if r > 0.01. 
On the contrary, if the cosmological B-mode signal is fainter, 
some strategy will be required to subtract the confusion 
noise produced by radio sources with flux-density S < 1 Jy. 
At sub-mm wavelengths, where radio sources become less 
and less relevant, confusion noise of point sources may be 
dominated by dusty galaxies (e.g., [79 1). Polarization of these 
objects should be relatively low, but they are expected to be 
significantly clustered. Although the level of this effect is still 
quite uncertain, polarization spectra from dusty galaxies could 
begin to dominate the one produced by ERS already at ^>200- 
300 GHz, at least at angular scales relevant for the detection 
of the B-modes. 
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